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Using time-domain Terahertz spectroscopy we performed direct studies of the photoinduced sup- 
pression and recovery of the superconducting gap in a conventional BCS superconductor NbN. Both 
processes are found to be strongly temperature and excitation density dependent. The analysis of 
the data with the established phenomenological Rot hwarf- Taylor model enabled us to determine 
the bare quasiparticle recombination rate, the Cooper pair-breaking rate and the electron-phonon 
coupling constant, A = 1.1 =b 0.1, which is in excellent agreement with theoretical estimates. 
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Soon after the first tunneling experiments in super- 
conductors (SCs) revealed that the quasiparticle (QP) 
tunneling can be qualitatively described in terms of the 
QP energy band picture, in which the SC can be treated 
similarly to a narrow gap semiconductor, it was realized 
that SCs can be used as detectors for far-infrared light 
[l| . Since real- world applications of SCs require an under- 
standing of their properties under non-equilibrium con- 
ditions several attempts were made already in 1960's to 
determine the time scales and processes that govern the 
recovery of the SC state. It was soon realized [2] that 
the recombination process, where two QP recombine to 
form a Cooper pair, was dominated by the emission of 
phonons with hw > 2 A, A being the SC gap. However, as 
pointed out by Rothwarf and Taylor [31] , the re-absorbtion 
of huj > 2 A phonons leads to the so called phonon bottle- 
neck, where the recovery of the SC state is not governed 
by the bare recombination of two QPs into the conden- 
sate, but rather by the decay of hu > 2A phonons. The 
non-equilibrium state can be generated either by QP in- 
jection in tunnel junctions [H, or by excitation with pho- 
tons of energy larger than 2A. In particular, due to the 
rapid development in the generation of ultrashort laser 
pulses and related femtosecond pump-probe techniques, 
the latter approach has enabled direct studies of carrier 
relaxation dynamics with femtosecond time resolution. 
Most of the work reported to date focussed on high-Tc 
SCs; predominantly on cuprates (3-[l5| but more recently 
also on pnictides 16-19|. The main focus of the research 
was on the influence of the superconducting gap and the 
normal state pseudogap on the carrier relaxation dynam- 
ics [IB, 0, H ll^S One of the still controversial 
issues with cuprates is the effect of the d-wave gap on the 
relaxation phenomena, and the role of phonons, i.e. the 
or absence I 



existence 
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20| or absence |8|, |9|1 of the phonon bot- 
tleneck in this class of SCs. Recently, systematic studies 
of the photoinduced melting of SC have been performed 
in cuprates |lll-ll4|. Moreover, manipulation of the order 
parameter in Lai.84Sro.i6Cu04 by applying intense Tera- 



hertz (THz) electric fields along the c-axis has also been 
demonstrated [23]. Despite numerous studies in high- 
Tc SCs, whose ground state properties are not well un- 
derstood, only a few experiments on conventional BCS 
SCs exist |24l427j. but no systematic study of dynamics 
as a function of temperature (T) and the absorbed en- 
ergy density (^4). Studies of non-equilibrium dynamics in 
conventional SCs are an important first step to under- 
standing relaxation dynamics in more exotic SCs. Since 
for conventional BCS SC A lies in the low-THz range, 
and can be resonantly probed by time domain THz spec- 
troscopy (TDTS), the data interpretation should be less 
ambiguous. 

In this Letter we report on the first detailed study of 
the SC state relaxation phenomena in a conventional SC 
NbN [m over wide range of T and ^. Utilizing the TDTS 
we have studied the T-dependence of its complex conduc- 
tivity, cr(co'), as well as its T- and ^d-dependent dynamics 
following photoexcitation with a fs optical pulse. We 
show that (j(uj) can be well fit to the prediction of the 
BCS theory [29], enabling direct studies of the temporal 
evolution of A. Both, the Cooper pair breaking (CPB) 
and the SC state recovery were found to be T- and A- 
dependent, and could be well explained by the Rothwarf- 
Taylor (RT) phonon-bottleneck model [i,[2o|. From the 
dependence of the CPB on A we were able to determine 
the microscopic QP recombination rate, i?, as well as the 
value of the electron-phonon (e-ph) coupling constant, A. 

NbN thin films were deposited by dc magnetron sput- 
tering on MgO substrates ^]. Films with thicknesses 
between 10-15 nm and T^ between 14.3-15.4 K (AT^ = 
0.19 - 0.17 K) were investigated. The TDTS set-up 
was built around a 250 kHz amplified Tiisapphire laser 
system, utilizing large area interdigitated photoconduc- 
tive emitter for the generation of THz pulses |32| , while 
the THz electric fields transmitted through the sample, 
Etr{t')^ and reference, Ere{t')^ were detected using the 
Pockels effect in GaP. The real, cri(cj), and the imagi- 
nary, (12(0;), parts of the optical conductivity were deter- 
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mined using the appropriate Fresnel equations. In non- 
equilibrium experiments, the films were excited by 50 fs 
pulses at the carrier wavelength of 800 nm. The time evo- 
lution of the SC state as a function of td^ the time delay 
of the THz pulse with respect to the optical excitation, 
was studied by either directly measuring aiuj^td) or by 
measuring the induced changes in the transmitted elec- 
tric field, AEtr{t\ td), at a fixed t'. The latter, spectrally 
integrated studies, are particularly useful for studying 
dynamics at low excitations 0, [l, [Isl • The THz beam 
diameter on the sample was ^1.5 mmjsol], while that of 
the photoexcitation beam was 4 mm, ensuring the homo- 
geneous lateral excitation profile. From reflectivity and 
transmission of the films at 800 nm we determined the 
optical penetration depth I opt ^12 nm. The absorbed 
energy densities at the film surface, ^4, were calculated 
from the reflectivity and lopt and are given in mJ/cm^. 



is at each T normahzed to its peak value (in a BCS SC, 
the transmissivity ratio peaks just above 2 A fsS^). 
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FIG. 1: (color online) The T-dependence of the a) imaginary 
and b) real part of cr (uu) of 15 nm NbN film on MgO sub- 
strate (Tc = 15.4 K). Solid lines are fits with the BCS equa- 
tions [29]. Inset to a) shows the transmitted THz transients 
through the NbN film below and above Tc; the arrow denotes 
the time t' = to with maximum change in the transmitted 
electric field. Inset to b): the T-dependence of A, extracted 
from fits to a (u) (symbols) overlaying the normalized trans- 
missivity ratio, TR. The BCS T-dependence of A is shown by 
the dashed line. 



Fig. 1 presents the T-dependence of a{uj) of a 15 
nm NbN film with Tc = 15.4 K. In the normal state 
a{uj) is fit with the Drude model with the plasma fre- 
quency Up = 460 THz (15300 cm~^) and the scattering 
rate h/r = 264 cm~^, in good agreement with studies 
on thick films [33|. In the SC state an excellent agree- 
ment between the data and a{uj) using the BCS model 
with a finite normal state scattering rate [29] (solid lines) 
is obtained. The extracted T-dependence of A, with 
A(0) ^ 0.75 THz (3.07 meV) and 2A/kBTc ^ 4.6, is 
shown in inset to panel b). A(T) is plotted on top of 
TR, the transmissivity ratio T(cj, T)/T(a;, 16K), which 
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FIG. 2: (color online) a) The time evolution of normal- 
ized transmissivity ratio (TR), recorded at A = 22 mJ/cm^, 
together with A (t) extracted from ai{(jO,td) fit with the 
BCS formula (2^1. b) Comparison of A(td) extracted from 
o'i{oo,td) (symbols) and from the spectrally integrated re- 
sponse AEtr{to,td) (dashed line). Inset: the measured T- 
dependence of AEtr{to) used to relate AEtr{to,td) to the 
change in A. 



Fig. 2a) presents the temporal evolution of TR follow- 
ing photoexcitation with a 50 fs optical pulse with A = 
22 mJ/cm^. The data reveal a strong suppression of SC 
on the 10 ps time scale followed by the recovery on the 
100 ps timescale. On the same plot A (t^), obtained by 
best fit to ai{uj^td) using the BCS model [29| is plotted 
by solid symbols. 

To perform systematic studies of the time evolution of 
A as a function of A and T, avoiding the signal drifts as- 
sociated with the long term laser stability, we studied the 
dynamics of the induced changes in the transmitted elec- 
tric field, AEtr(t'), at a fixed point of Etr{t' = to). The 
fixed point of Etrito) in these experiments was at the 
point of maximum time derivative of the electric field, 
to = ps, where the changes in Etr{to) upon entering 
the SC state are the largest — see inset to Fig. la). In 
order to obtain a direct link between AEtr{to) and the 
photoinduced change of the gap, (5 A, we first measured 
AEtr{to,T, To) = Etr{to, T) - Etr{to, To), where To is the 
T of the sample at which photoinduced studies are per- 
formed. Combining A£^^r(to, T, 4.3K), shown in inset to 
Fig. 2b), with the T-dependence of the gap (inset to 
Fig. lb)) enabled us to extract the temporal evolution 
of the gap from AEtr{to,td) traces. The trace obtained 
this way at A = 22 mJ/cm^, shown by the dashed curve 
in Fig. 2b), matches well the time evolution of the gap 
extracted from the ai{uj,td) data. Therefore, by mea- 
suring the spectrally integrated response, the T- and A- 
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dependence of the photoinduced gap change, ^A, can be 
extracted quickly, avoiding system drifts. 
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FIG. 3: (color online) a) The relative change in gap, ^A/Ao 
recorded at 4.3 K for various A in mJ/cm^. The dashed lines 
are fits to the data with Eq.(l). b) The dependence of ^A/Ao 
on at 4.3 K. The dashed line is a fit to the simple saturation 
model, the solid line is the linear fit. Inset: The T-dependence 
of the saturation energy density, Asat^ compared to the T- 
dependence of A^ (dashed line). 



Fig. 3a) presents the SA{td) traces for various A 
recorded at 4.3 K. As shown in Fig. 3b) the maxi- 
mum induced change initially increases linearly with ex- 
citation density, followed by a saturation resulting from 
suppression of SC. To estimate the characteristic A re- 
quired to suppress SC, we use a simple saturation model, 
SA/Aq = {1 — ex.-p {—A/Asat))- Here Agat is the char- 
acteristic absorbed energy density required to suppress 
SC and is Asat (4.3 K) ^ 25 mJ/cm^. We find that 
in NbN Agat is, within the experimental accuracy, com- 
parable to the thermodynamic SC condensation energy, 
Ec = 22 mJ/cm^. Here Ec = = N{0)A'^/2, with 

Be = 0.234 T [m being the thermodynamic critical field 
and N{0) the single spin density of states at the Fermi 
level. As shown in inset to Fig. 3b), the T-dependence of 
Asat is found to follow the T-dependence of A^ (dashed 
line). These observations are in strong contrast to sim- 



ilar studies in cuprates ll,[l3|, [l^, where Asat is found 
to be about one order of magnitude higher than Ec and 
nearly T-independent Q. 

As demonstrated in Fig. 3a), both the CPB and the SC 
recovery show pronounced dependence on A. The CPB 
rate is found to increase with increasing A^ similarly to 
the case of MgB2 [26]. This behavior can be attributed 
to the intrinsic non-linearity of the SC state relaxation 
process, where the populations of the photoexcited QP 
and the high frequency {hw > 2A) phonon (HFP) den- 
sities are described by the coupled RT rate equations 
[1, [20I, [Isl . Here the long timescale of the CPB implies 
that photoexcited hot electrons (holes) initially generate 
a high density of HFP, which subsequently break Cooper 



pairs until the q uasi-e q uil ibr ium between QP and HFP 
populations is reached [20|, |26| . 

The CPB can be quantitatively analyzed using the 
RT model. In the low photoexcitation limit {6 A << 
A), where the RT analysis is applicable, ^A is propor- 
tional to the photoexcited QP density, npj. It is easy 
to show that at 4K npj is for the entire range of A 
substantially higher than the density of thermally ex- 
cited QP, riT. From ut ^ A^(0)v'2^rM^exp(-A/T), 
where 7V(0) = 0.44 spin-^unit ceir^ eV-\[36] it fol- 
lows that riT (4 K) ^ 1 10~^ unit cell"^ while A = 
0.14 mJ/cm^ corresponds to npi ~ 3 10~^ unit cell~^. 
Therefore, the CPB {td < 20 ps), which is well separated 
from the SC state recovery dynamics (^100 ps), can be 
fit (dashed lines in Fig. 3a)) with ^ 



npi (td) 



R 



4 2r rl-Kexp{-tdP/r) 



. (1) 



Here /3 is the CPB probability by absorption of HFP, 
R is the bare QP recombination rate, while K and r 
are dimensionless parameters (Eq.(3) of Ref . (26| ) deter- 
mined by /3, R and the initial conditions (the ratio of 
the absorbed energy in the HFP and QP subsystems fol- 
lowing the initial e-e and e-ph scattering of hot electrons 
(holes)), which should be independent on A for these low 
excitation densities 2o|, [26| . The extracted dependences 
of r//3 and i^T on ^4 are shown in Figs^ 4a) and 4b), to- 
gether with the fit with Eqs.(3) of Ref. [26] (dashed lines). 
The best fit is obtained when 91 % of ^4 is initially trans- 
ferred to the HFP subsystem, giving the values of the 
microscopic constants = 6 ± 1 ps, = 160 ± 20 
ps~^unit cell. 

Fig. 4c) presents the dependence of the SC recovery 
time, Trec^ on T and A. Here Tree is obtained by fitting 
the recovery dynamics with a single exponential decay. 
As shown in inset to Fig. 4c), at 4 K r^J. first increases 
linearly with A^ mimicking the intrinsic bimolecular ki- 
netics of the QP recombination [20]. For ^ ^ 5 mJ/cm^ 
and for T > 10 K r^J. is constant at t~^].{0) ~ 0.01 ps~^. 
In NbN T^c(O) i^ governed by the escape of the HFP 
into the substrate and is inversely proportional to the 
film thickness [s^ . It was argued that the dependence of 
Tree on ^7 obscrvcd in cuprates, implies the absence of 
the phonon bottleneck in cuprates [H, ^ . The subsequent 
analysis of the RT model [20[, however, suggested a (per- 
haps counterintuitive) result, that also in a phonon bot- 
tleneck case r^J. oc npi for < npi < /3/R. Here /3/R 
is the material dependent characteristic QP density, be- 
low which the RT perturbative description is applicable 
(20). Indeed, A = b mJ/cm^ corresponds to np/ ~ 0.001 
unit cell~^, identical to P/R ~ 0.001 unit cell~^ deter- 
mined from the analysis of the CPB. Similarly, the fact 
that r~^l is above ^ 10 K A- and T-independent follows 
from the fact that nT(lOi^) ~ 0.0003 unit cell~^ is com- 
parable to P/R. Therefore, as experimentally demon- 
strated here for a standard BCS SC, the 74-dependent 
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FIG. 4: (color online) Panels a) and b) show the dependence 
of r//3 and K on A; the values are extracted from fits to the 
Cooper pair breaking dynamics (Fig. 3a). Panel c) shows the 
T-dependence of the SC recovery time, Tree, for several A (in 
mJ/cm^). Inset: T~el {^) recorded at 4.3 K. For low A the 
relaxation rate increases linear with A (dashed line). 



SC state recovery at low Ts seems to be an intrinsic fea- 
ture of a SC driven out of equilibrium, and should be 
observed in the phonon bottleneck case, providing that 
riT < Up I < P/R. Based on this study, we argue that 
it is the absence of ^d-dependence of the SC state recov- 
ery, observed e.g. in optimally doped YBCO [H^ and in 
an over doped BSCCO [8], that presents an anomalous 
behavior and not vice versa as argued (8|, l9| . 

As demonstrated, the dynamics of photoexcited NbN 
can be over large range of T and A well described by the 
phenomenological RT model, enabling us to determine 
the values of microscopic parameters R and p. Since R = 
hN (o)Q'^ im , where Qc is the phonon cut-off frequency, 
the value of A can be determined. Taking the known 
values for A, 7V(0) and = 16 THz [28J, we obtain 
A = 1.1±0.12, which is in a very good agreement with the 
theoretical estimates, A = 1 — 1.12 36|, 37|. Importantly, 
such an approach does not suffer from the (ambiguous) 
underlying assumptions of the two-temperature model, 
the model commonly lo|, [Hj used to determine A. 

In summary, we presented the first systematic studies 
of the SC state relaxation phenomena in a conventional 
BCS superconductor NbN with ps time-resolution. Uti- 
lizing the TDTS, we were able to study the time evo- 
lution of A over large range of temperatures and excita- 
tion densities. We demonstrated, that both the CPB and 
SC state recovery dynamics depend strongly on the ex- 
citation density in agreement with the predictions of the 
phenomenological RT model [^,|26|. Studying the CPB 
enabled us to determine the values of the microscopic pa- 
rameters i^, /3, as well as the dimensionless e-ph coupling 



constant A in NbN. This approach could be used to de- 
termine the electron-boson coupling strengths in high-Tc 
SCs, providing the CPB dynamics as a function of exci- 
tation density can be experimentally resolved. Last but 
not least, our results on NbN emphasize the unusual na- 
ture of the 74-independent relaxation observed in some 
cuprates [I, ^ . 
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